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ABSTRACT 

Black holes exceeding a billion solar masses have been detected at redshifts greater than six. The 
rapid formation of these objects may suggest a massive early seed or a period of growth faster than 
Eddington. Here we suggest a new mechanism along these lines. We propose that in the process of 
hierarchical structure assembly, dense star clusters can be contracted on dynamical time scales due 
to the nearly free-fall inflow of self-gravitating gas with a mass comparable to or larger than that of 
the clusters. This process increases the velocity dispersion to the point where the few remaining hard 
binaries can no longer effectively heat the cluster, and the cluster goes into a period of homologous core 
collapse. The cluster core can then reach a central density high enough for fast mergers of stellar-mass 
black holes and hence the rapid production of a black hole seed that could be 10 5 M Q or larger. 
Subject headings: galaxies: formation — galaxies: evolution — galaxies: nuclei — black hole physics 



1. INTRODUCTION 

Evidence from observations of high redshift ac- 
tive galactic nuclei (AGN) shows that massive black 
holes (m bh ~ 10 9 M^) have alread y formed at 
red shifts z > 6 (jWillott et al.l 120031 : iBarth et al.l 
2003). It is challenging to produce such masses so 
early via Eddington-limited gas accretion onto ini- 
tially stellar-mass black holes, hence there have been 
various explorations of the possi bility of higher-mass 
seeds, e.g., Populati on III stars (jMadau fc Reesl 120011 : 
iVolonteri et al.ll2003lL massive stars formed through run- 
away collisions dDevecchi fe Volonteril I2009D or quas i- 
stars (jBegelman et al.ll2006t lLodato fc Nataraianl l2006). 
Here we renew a earlier suggestion: that a sufficiently 
massive and dense star cluster with a central black 
hole sub-cluster could undergo collapse, including merg- 
ers due to gravitational radiation, building up a mas- 
sive black hole much faster than would be possible by 
Eddington-limited gas accretion onto a stellar-mass seed. 

SMBH formation within stellar clusters in galactic 
nuclei has been explored in a ser ies of papers by 
iQuinlan fc Shapiro! (fl987l . [l989l [l990h . They considered 
core collapse in a cluster of compact objects where the 
cluster core contracts driven by two-body relaxation to 
the point where the central potential well is so deep that 
a relativistic instability sets in and a single black hole 
is formed via the collapse of the central core of compact 
objects. However, this instability requires a central red- 
shift z ~ 0.5, meaning that the core would have to be 
extremely massive and dense. 

In this Letter we consider a variation on the cluster- 
collapse model. We begin with cluster conditions simi- 
lar to those seen in nuclear stellar clusters, which have 
masses a nd radii roughly comparable to globular clus- 
ters (e.g. JSeth et al.ll2008t ). In these clusters, the more- 
massive stellar-mass black holes are likely to segregate 
to the cluster center forming their own dark core. Bi- 
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naries within the core will provide the cluster with en- 
ergy via binary-single encounters. These binaries are rel- 
atively wide (~ 0.1 — 1 AU), hence the timescale for them 
to merge via the effects of gravitational radiation is ex- 
tremely long. Binary heating will thus provide a fuel to 
support the cluster core from complete collapse. 

Infall of gas into such a nuclear stellar cluster may 
have profound consequences for its evolution. Infall is 
likely to occur du ring a merger between two galaxies 
dMaver et al.ll2010f). which w ill be common at high red- 
shifts ( Bellovary et al.l 120111 ) . The addition of signifi- 
cant mass will cause the cluster to shrink after virial- 
ization, and in some circumstances significant accretion 
onto cluster stars may occur. The cluster potential well 
will be deepened, increasing velocity dispersions and de- 
creasing the number and semi-major axes of the binaries 
that heat the cluster. Because the inspiral timescale due 
to gravitational radiation depends sensitively on the ini- 
tial separation, these smaller binaries will have signifi- 
cantly decreased merger timcscalcs. Mergers remove bi- 
naries and thus remove the energy source for the cluster. 
Without an energy source, the cluster core will undergo 
deep core collapse, potentially causing a runaway merger 
of the black holes residing there. For this collapse to 
occur we require not only that the merger timescale of 
binaries due to the emission of gravitational radiation 
is small, but also that the merger products (which will 
receive kicks due to the asymmetric emission of gravita- 
tional waves) are retained within the cluster core. 

In Section 2 we set out the conditions for a nuclear 
stellar cluster such that black-hole binaries within the 
core are likely to merge (and be retained) on shorter 
timescales than cluster heating via binary-single encoun- 
ters. In Section 3 we consider the immediate effects of in- 
fall gas into a nuclear stellar cluster. In Section 4 we con- 
sider the subsequent evolution of a nuclear stellar cluster 
after an episode of gas infall, possibly leading to a chain 
of events producing an SMBH. 

2. DRY NUCLEAR STELLAR CLUSTERS 

We begin by considering a cluster free of gas and as- 
sume the cluster contains only two species: stellar-mass 
black holes, and less-massive main-sequence stars. The 
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heavier black holes will sink to the core by the effects 
of mass segregation. The cluster will likely be vulner- 
able to the Spitzer instability, whereby the black holes 
will form their own central sub-cluster within the cluster 
core. Although some black hole - black hole interactions 
can lead to ejections, a significant fraction of the initial 
black hole population rem ains even after several Gy r at 
standard cluster densities (jMackev et al.ll2007l[2008[) . 

Black hole binaries within the dark core that are soft 
(ie their binding energy is less than the typical kinetic en- 
ergy of a single black hole) will be split rapidly by binary- 
single interactions. In contrast, hard binaries tend to 
harden further after binary-single interactions (and heav- 
ier black holes tend to swap into the binaries), meaning 
that the single and binary both get a kick after the inter - 
action and the cluster is heated as a result (Heggie ll975[ ). 
The timescale for a given object to have a gravitationally 
focused enc ounter with another object within the cluster 
is given by (jBinnev & Tremaindl2008l ) 

w = 7 x lQ w nz 1 Voc,wr~ t m~ 1 yr (1) 

where n§ is the number density of stars/black holes in 
units of 10 5 objects/pc 3 , r mm is the minimum distance 
during the encounter in solar radii, m is the sum of the 
masses of the objects involved in the encounter in solar 
masses, and i>oo,io is the relative speed at infinity in units 
of lOkm/s. When considering encounters between bina- 
ries and single stars, it is reasonable to set r m - m = a\,- m , 
the semi-major axis of the binary. If we consider en- 
counters between binary and single black holes, all of 
equal mass, then m = 3mt,h- The number density of 
objects in the core is n ~ 3-/Vbh/47rr- 3 , where Abh is the 
total number of black holes and r c is the core radius. 
Assuming the cluster is in virial equilibrium, one may 
estim ate that Vqc ~ y/QAGM c /rh (jBinnev fc Tremamel 
2008), where M c is the total cluster mass. This implies 

Woo, 10 — 4.36a/Mc fi/ r h where M c $ is the total cluster 
mass in units of 10 6 M Q and is the cluster half-mass 
radius in pc. 
The semi-major axis of a binary is given by 
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where x is the ratio of the binary binding energy to the 
kinetic energy of the stars/black holes, and mbh is given 
in solar masses. One can show that the timescale for an 
encounter between a given black-hole binary and a single 
black hole is given by 
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where r c is the cluster radius in pc. Setting r c = fir^, 
nibh = IOM0, m* = IMq, and Abh = aJV t , where A* is 
the total number of stars in the cluster, we obtain 
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Assuming a Salpeter IMF, with minimum stellar mass 
of 0.2 M Q and that all stars more massive than 25 M Q 



produce black holes, we have a 
gives 
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Tight binaries will lose energy via the emission of grav- 
itational radiation. The timescale for two black holes, in 
a binary of initial separati on abin and e ccentricity e, to 
spiral together is given by (|Peterslll964h 
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Binary-single scattering will leave the binaries with a 
thermal distribution of eccentricities where the distribu- 
tion of orbital eccentricities is given by dn/de — 2e. The 
median eccentricity is thus e mc d = 1/V2, hence a typical 
binary merger time is reduced by a factor ~ 10 and could 
be reduced by much more for higher eccentricities. Using 
the expression for abin given above, we obtain 
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Fig. 1. — Log of timescales (in years) for binary-single encoun- 
ters (dashed line) and gravitational radiation inspiral (solid line) 
in binaries as a function of log of velocity dispersion within clus- 
ters, t>oo lOi in units of lOkm/s, assuming binaries with hardness 
parameter x = 3, and eccentricity e = l/v2. The vertical dotted 
line represents the value of i>oo,lo when merging black holes are 
likely to be retained. 



We plot r gr and T2+1 as a function of Woo.io in Fig- 
ure [l] taking x — 3, and assuming the binaries initially 
have the median eccentricity of a thermal distribution 

(e = e mc d = 1/V2). Thus, for moderately hard (x ~ 3) 
binaries in a typical nuclear stellar cluster (u^io ~ 3), 
i gr is extremely large: black hole binaries are unlikely 
to merge via the effects of gravitational radiation. If 
r gr < T2+1, binaries will merge before they have the 
chance to heat the cluster. The exact value of i>oo,io 
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when this is the case will depend on the hardness and 
eccentricity of the binary, but for moderately hard, and 
eccentric, binaries, this is likely to be true for Doo,io r& 50. 

When two inspiraling black holes merge, they receive 
a kick due to the asymmetric emission of gravitational 
radiation. This kick depends on the mass ratio of the 
two black holes and their spins, but can be as large 
as ~ 4000 km s" 1 for an optimal configuration (e.g., 
iBaker et al.l l2008h . In order to retain merger products 
within the stellar cluster, we require the cluster escape 
speed to exceed 1000 km/s (or Uoo,io ~ 50). This is 
equivalent to requiring that 
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Any merger products ejected extremely close to the es- 
cape speed may be left on very wide orbits outside of 
the cluster, however objects ejected at less than about 
80 % of the escape speed will remain in the cluster. We 
thus see that the conditions required for the retention of 
merger products are rather similar to those which give 
us T gr < T2+i- If the merger products are retained within 
the cluster, we expect the dark core to continue to con- 
tract, having the potential to lead to a runaway merger 
of black holes. The outcome of such a process could be 
the production of an intermediate-mass black hole which 
would act as the seed for an SMBH. 

However, most observed nuclear stellar clusters will not 
satisfy equation (|SJ|. We therefore require some mecha- 
nism to cause the cluster to shrink and/or increase its 
mass in order to reach the conditions required for run- 
away black hole mergers to occur within the dark core. 

3. THE EFFECT OF GAS INFALL INTO THE CLUSTER 

Suppose that gas with mass comparable to or larger 
than the stellar mass of the cluster falls quickly to the 
center of the cluster. Just as in the inverse problem (e.g., 
where mass is expelled from a binary in a supernova) , we 
assume that the angular momentum of individual stellar 
orbits is conserved, meaning that the characteristic size 
of the orbits is proportional to 1/M, where M is the 
mass interior to the orbit. Thus, for example, adding a 
gas mass equal to the cluster mass halves the sizes of the 
orbits exterior to the gas and increases the stellar den- 
sity by a factor of eight. As shown below, at redshifts 
z > 10 one can indeed have gas inflows with total masses 
up to ~ 10 times the mass of the stellar cluster. The 
key question is whether much of the gas flows to the in- 
ner ~pc or closer, rather than stalling and forming stars 
much farther away, which would not contract the central 
cluster. 

One might imagine that a cooling flow could form 
that would bring the g as in. Indeed, s uch a scenario 
has been suggested by iVesperini et al.l (|201Ct) to grow 
intermediate-mass black holes in early globular cl usters . 
However, as demonstrated by iD'Ercole et al.l ([2008), 
the distributed injection of even a comparatively small 
amount of luminosity (Qcr ~ 6 x 10 37 erg s _1 in their 
treatment) is sufficient to hold off the cooling flow. This 
is considerably less than the Eddington luminosity from 
even a single 10 M Q black hole, and given that this lumi- 
nosity is typically distributed over a large spread in pho- 
ton energies that is therefore absorbed over a wide range 



in gas column depths, it appears unlikely that a cooling 
flow will develop in the situation we consider here. 

We instead tu rn to the scenario proposed by 
IMaver et al.l ((2010), in which the self-gravitating gas is 
subject to instabilities that funnel much of the gas to 
the center in a low angular momentum flow that, in 
their calculations, gets to 0.2 pc or closer to the cen- 
ter. Gravitationally driven turbulence prevents the gas 
from fragmenting into stars. Such a flow may lead to 
the fo rmation of a massive black hole, as IMaver et al.l 
(2010) suggest, but the large range of scales bet ween the 
final ~ 0.1 pc simulated by IMaver et al.l [20101 and the 
~ 0.01 AU scale of the black hole means that the final 
collapse is still unresolved and might be hampered by a 
variety of effects. Nonetheless, the inflow will also very 
effectively contract an existing stellar cluster, leading to 
high central densities, rapid mass segregation, and fast 
interactions among the existing stellar-mass black holes 
that could lead to quick coalescence and the formation 
of a massive black hole seed. This scenario is particu- 
larly probable given that such processes likely occurred 
in miniature in the smaller clusters and gas flows that 
accumulated to form the central cluster. Thus the ini- 
tial black hole masses would have had a wide range; this 
is known to speed up ma ss segregation and coalescence 
(jQuinlan fc Shapiro! Il987l ). 
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Fig. 2. — Amount of mass as a function of time for gas (black 
line) and stars (dashed line) in the cosmological simulation. Masses 
are the total mass enclosed within 520 pc (two softening lengths) 
from the center of the primary galaxy in the simulation. 



As an illustration of the gas infall into nuclear regions 
which can occur, we turn to high-resolution cosmological 
simulations of galaxy formation, considering the evolu- 
tion of a massive galaxy destined to become a massive 
elliptical at low redshift (|Bellovarv et al.l [201 If ). In Fig- 
ure [2] we show the mass enclosed within 520 pc of the 
galaxy center vs time for the gas (black line) and stars 
(dashed line). The gaseous and stellar inflow are due to 
a combination of direct accretion of matter from the cos- 
mic web/filaments as well as galaxy mergers. From this 
plot we see that significant inflow of both gas and stars 
occurs at high redshift. Repeating this measurement for 
similar simulations reveals that high redshift (z > 10) 
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massive inflow on this scale is a common occurrence for 
galaxies in this mass range (a few xlO 11 M Q at z = 5). 

One might wonder whether the rapid contraction of a 
cluster by a factor of a few would lead to a high rate of in- 
teractions of single objects with hard binaries and hence 
to re-expansion of the cluster due to energy input. The 
degree to which this can happen clearly depends on the 
fraction of stars or black holes in hard binaries that are 
not so hard that they merge or collide rapidly. This frac- 
tion, in turn, depends on the velocity dispersion; higher 
velocity dispersion means fewer hard binaries and less 
binary binding energy that can potentially be tapped to 
hold off core collapse. We also note that close three-body 
interactions between objects of comparable mass yield a 
thermal distribution of eccentricities. The eccentricities 
that give pericenter distances low enough for collisions 
(~ 0.01 AU for solar-type stars) or fast merger by grav- 
itational radiation (also ~ 0.01 AU for ~ 10 Mq black 
holes to merge in a few million years or less) will destroy 
the binary. If the hard-soft boundary is a < 1 AU (corre- 
sponding to a velocity dispersion a ~ 30 km s _1 for solar- 
type stars or a ~ 100 km s _1 for ~ 10 Mq black holes), 
the available binding energy in binaries is less than the 
binding energy in the singles, hence binary-single inter- 
actions are inefficient at holding off core collapse. 

4. THE FATE OF THE MASS-LOADED CLUSTER 

We now consider the subsequent evolution of the stellar 
cluster, assuming it has received a significant influx of 
gas, and perhaps stars. Beginning with a nuclear stellar 
cluster of mass M C) 6 ~ 1, and half-mass radius ~ 1, 
the infall of 10 7 M© of gas is likely to shrink the cluster 
by a factor of about ten leaving initial core densities n ~ 
10 s stars pc -3 and velocity dispersions «oo,io ~ 30 — 50. 
Several things will happen to this cluster. Wider binaries 
(which previously heated the cluster and prevented core 
collapse) will now be soft, as the velocity dispersion has 
increased, and will be broken up. Any remaining binaries 
will quickly merge via gravitational radiation inspiral or 
collisions. Merger products are likely to be retained with 
the cluster owing to the increased escape speed which 
resulted from the gas infall. It is possible that some 
merger products will in turn exchange into other binaries 
which then merge. This process may even be repeated. 

After a short time, the cluster core will contain black 
holes having a broader range of masses than previously. 
In addition, all binaries will either have been broken up 
or have merged. Thus the cluster core will lose its source 
of energy, and will begin to contract approaching core 
collapse as it transfers energy to the cluster halo via 
two-body scattering. As the core density increases, run- 
away mergers between any no n-compact stars within th e 
cluster core becomes possible (jQuinlan fe ShapiroH l990). 
Such collisions may produce more massive stars which 
in turn will add to the population of black holes and 
neutron stars. 

At high densities, extremely close encounters occur be- 
tween two compact objects where binary formation via 
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gravitational wave emi ssion is possible. The tim escale 
for capture is given by (|Quinlan fc Shapiro! 119891) 

r pb 7 x lO 12 ^ 1 /T 2/7 'm~ 12/ V 11/7 vr (9) 

leap ~ I A 1U H 5 fl w JM 10 o y oo,500 y 1 ' V y J 

where /iio is the reduced mass of the two compact objects 
capturing each other (in units of 10 M Q and Afioo is their 
total mass (in units of 100 Mq). We see that when the 
core reaches densities n ~ 10 12 stars pc -3 , which at a ve- 
locity dispersion v ~ 300 km s _1 contains few x 10 4 Mq, 
the formation of BH-BH binaries via gravitational wave 
emission becomes possible on a timescale <]Myr. These 
binaries will have very short inspiral timescales. We will 
thus have a second phase of black hole binary mergers. 
The black holes produced via such mergers have greater 
mass, and given that the timescale for mergers scales as 
M~ 12 ' 7 , higher -mass black holes have a greater probabil- 
ity of merger and the process will run away. Such a phase 
may well collect most of the mass of the stellar black holes 
into a single black hole, which would have several percent 
of the original stellar mass of the system for typical initial 
mass functions. Thus the mass of the single black hole 
could be ~ 10 5 M Q for situations similar to that shown 
in Figure [5] (assuming most infalling gas reaches the cen- 
tre). Accretion of neutron stars and white dwarfs, both 
of which will be swallowed whole for M ^ 10 5 Mq, could 
boost the mass by a factor of a few, depending on the 
mass at which stellar interactions in the radius of influ- 
ence of the black hol e become an effe ctive heating source 
for the cluster (e.g.. IGill et all 120081) . Overall, the time 
elapsed between mass infall to the pr oduction of the seed 
SMB H will be roughly 100 Myr fe.g.. lQuinlan fc Shapiro! 
1990), hence producing a ~ 10 5 Mq black hole within 
~ 300 — 400 Myr after the beginning of the universe. 
Growth to ~ 10 9 Mq via gas accretion could then oc- 
cur by a redshift z ~ 6 via standard Eddington-limited 
accretion onto moderately spinning black holes. 

In summary, we have proposed that clusters with initial 
structural parameters similar to current-day nuclear clus- 
ters may (1) form as part of early hierarchical merging, 
(2) accrete gas with a total mass comparable to or greater 
than that of the cluster at redshifts z > 10, (3) contract 
as a result so that binary heating is ineffective, (4) un- 
dergo core collapse to a density high enough that stellar- 
mass black holes merge, and thus (5) have most of the 
mass originally in stellar-mass black holes collect in to a 
single black hole that could be ~ 10 5 Mq or larger. This 
black hole would therefore be a high-mass seed that could 
comfortably grow to supermassive size by the observed 
redshifts z ~ 6. 
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